Cooling many particles at once 



Almut Beige/'^ Peter L. Knight,^ and Giuseppe Vitiello^ 

^Department of Applied Mathematics and Theoretical Physics, 
University of Cambridge, Wilherforce Road, Cambridge CBS OWA, UK 
^Blackett Laboratory, Imperial College, Prince Consort Road, London SW7 2BW, UK 
^ Dipartimento di Fisica "E. R. Caianiello," LN.F.N. and LN.F.M., 
Universitd di Salerno, 84 100 Salerno, Italy 

February 9, 2008 

Abstract 

We propose a mechanism for the cohective cooHng of a large number A'^ of trapped 
particles to very low temperatures by applying red-detuned laser fields and coupling 
them to the quantized field inside an optical resonator. The dynamics is described by 
what appears to be rate equations, but where some of the major quantities are coher- 
ences and not populations. The cooperative behavior of the system provides cooling 
rates of the same order of magnitude as the cavity decay rate k. This constitutes a 
significant speed-up compared to other cooling mechanisms since k can, in principle, 
be as large as \/iV times the single-particle cavity or laser coupling constant. 

Pacs numbers: 03. 67. -a, 42.50.Lc 

1 Introduction 

Cooling and trapping techniques have improved dramatically over the last few decades. A 
very efficient method to transfer, for example, a single atom rapidly to a very low temperature 
is sideband cooling This requires a red-detuned laser field whose detuning equals the 
frequency of the vibrational mode of the atom. When the laser excites the system from the 
ground to an excited state, the vibrational energy reduces by one phonon. Afterwards, the 
atom most likely returns into the ground state via spontaneous emission of a photon and 
without regaining energy in the vibrational mode. The corresponding non-unitary evolution 
effectively reduces the temperature of the atom and yields an overall decrease of the von 
Neumann entropy in the system [2]. Other experiments aim at cooling molecules, which 
have a much richer inner level structure and are therefore harder to control than atoms. 

Here we apply the idea of side band cooling to a large number of particles (atoms, ions or 
molecules). As in the one-atom case, a red-detuned laser field excites the particles, thereby 
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Figure 1: Atomic level scheme (a) and experimental setup for the collective cooling of many 
particles trapped inside an optical cavity (b) or an optical ring resonator (c). 



continuously reducing the number of phonons in the system. To return the particles into 
their ground state, they should couple to the quantized field of a leaky optical cavity, whose 
frequency a;cav equals the dipole transition frequency ujq of each particle (see Figure HJ^a)) . 
Once the particles transferred their excitation into the resonator mode, it leaks out through 
the cavity mirrors. Using cavity decay, instead of spontaneous emission from excited levels, 
helps to avoid heating due to rescattering of photons within the sample. It also minimises 
spontaneous emission into unwanted states and allows to control even complicated level 
structures, like molecules. 

Crucial for obtaining maximum cooling is the generation of cooperative behavior of the N 
particles in the excitation step as well as in the de-excitation step. This is possible when the 
Rabi frequency of the laser field for the cooling of a vibrational mode with frequency u 
is for all particles the same and all particles see the same cavity coupling g (small variations 
of g and Qi, and non-ideal initial conditions do not substantially affect our conclusions and 
will be considered elsewhere [3]). Realising this inside an optical resonator requires self- 
organization of the particles in the antinodes of the cavity field, as predicted in Ref. [HE]. 
If the dipole moments of the atoms are in the average parallel to the cavity mirrors, the 
Rabi frequencies fly are practically and to a good approximation for all particles the same, 
if the laser field enters the cavity from the side as shown in Figure ^b). The laser could 
also enter the cavity through one of the cavity mirrors, as it was the case in the many-atom 
cavity QED experiment in Grangier's group in 1997 Alternatively, a ring resonator can 
be employed, as in Refs. [H El El EI , if a laser field with cjiaser enters the resonator in a certain 
angle^ ip with cos ip = o^iaser/i^cav (see Figure [Tfc)). 

We also require that the particles are initially all prepared in their ground state (in the 
large N limit fluctuations can be neglected). Then the time evolution of the system remains 
restricted within a Dicke-symmetric subspace of collective states fUj • As shown below, these 
states experience a very strong coupling to the laser field as well as to the cavity mode and 
the system evolves into a stationary state with no phonons on the time scale given by the 
cavity photon life time. Several schemes for the cooling of atomic ensembles have been 
proposed [TTl [T21 [121 HI IHj and first cavity-cooling experiments have already been performed 

^This assures that the phase factors in the cavity interaction term of the Hamiltonian can be absorbed 
into the definition of the excited states |f) of the particles in the respective positions. The same applies for 
the phase factors of the laser amplitudes, which are now exactly the same as the phase factors of the cavity 
field term. 
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[HlITlin]. Compared to these, the scheme we propose provides a significant speed-up of the 
cooling process if operated in a regime with 

K ~ VNg, \VN7]VI > r with Vt= (1) 

V 

Here n denotes the cavity photon decay rate, F is the spontaneous decay rate of a particle 
in the excited state and 77 is the Lamb-Dicke parameter characterizing the steepness of the 
trap. For sufficiently large A^, condition (P) can be fulfilled even if 77 ^ 1 and the system 
is operated in the bad cavity limit with k ^ j > F. Moreover, we remark that Eq. (Q) 
describes a strong damping regime in which the system can only accumulate a small amount 
of population in the excited states of the particles. 

The proposed cooling scheme can be used to cool a large number of particles very effi- 
ciently. It is therefore an interesting question, whether the described setup might be used 
for the preparation of Bose-Einstein condensates. Currently, these experiments mainly use 
evaporative cooling ^Jji] which systematically removes those atoms with a relatively high 
temperature from the trap. Consequently, only a small percentage of the initially trapped 
atoms is finally included in the condensate. If one could instead cool all the atoms efficiently, 
yet at the same time avoid the loss of particles, it should become easier to experiment with 
large condensates. Besides, cooling is also crucial for ion trap quantum computing where 
the achievable gate operation times can depend primarily on the efficiency of the cooling of 
a common vibrational mode ^H] • 



2 Bosonic behavior of a large atomic sample 

We consider a collection of two-level particles with ground states |0)j and excited states 
each of them described by a^i = — |0)ij(0|) with eigenvalues ±|. Transitions 

between the two levels are generated by = |1)m(0| and a~ = |0)jj(l|. Our fermion-like 
iV-body system is thus described by the su(2) algebra 

[a3,cr±] =±a±, [a-,a^] = -2as (2) 

with cr^ = erf and (T3 = Y^iLi'^^i- Under the action of o"^, describing the laser ex- 

citation, the initial state with all particles in the ground state, |0)p, is driven into the 
Dicke-symmetric states |/)p with 

|/)p = [\0,02...0n-i1n~i+i1n-i+2...1n) + ... + \hh...hOi+iOi+2...0N)]/^l{^ , (3) 

a superposition of all states with / particles in The difference between excited and 
unexcited particles is counted by since p(/|cr3|Z)p = I — jN. For any / 

CT+ |/)p = v/ITIv/iV^|/ + l)p, 

a-|/)p = y^iV-(/-l)v^|/-l)p, (4) 



3 



showing that and as are represented on |Z)p by the Holstein-Primakoff non-hnear boson 
reahzation pTj ITH] 

a+ = 4nS^As , a- = y/NAsS- , cxg = S+S' - ^A^ (5) 

with 

As = ^l-S+S-/N, 5+|/)p = Vm|/ + l)p, 5-|/)p = v^|/-l)p (6) 

for any /. The a's still satisfy the su(2) algebra (0). However, for N ^ I, Eqs. (@]) become 

a^\l), = VNS^\l), (7) 

and thus = a^/\fN for large N. In the large N limit, the su(2) algebra (0) written in 
terms of and 6*3 = contracts to the (projective) e(2) (or Heisenberg-Weyl) algebra 

[1211201 

[53,5±] = ±5±, [S-,S+] = l. (8) 

The meaning of Eqs. ^ and (|H1) is that, for large A^, the laser excites collective dipole 
waves, denoting the creation and annihilation operators of the associated quanta, and 
the collection of single two-level particles manifests itself as a bosonic system. 

3 Collective cooling of common vibrational modes 

Each particle may couple to its own phonon modes and there can also be common vibrational 
modes. We first discuss a scheme for the collective cooling of common modes. This requires 
the application of laser fields, each red-detuned by a phonon frequency v and with Rabi 
frequency VLy. In the following, is the annihilation operator for a phonon with v and c 
denotes the annihilation operator for a cavity photon. The Hamiltonian of the system in the 
interaction picture and within the rotating wave approximation^ then equals^'^ 

Hi = Y. l^^V^u S+K + hVNg S+c + H.c. (9) 

A detailed derivation of the atom-phonon coupling Hamiltonian in Eq. ^ can be found in 
Ref. It applies in the Lamb-Dicke limit, where the atom-phonon coupling is relatively 
small compared to the phonon frequency z/ and (^T^fi,^)^ <^ u"^. Moreover, we neglect the 

^This approximation introduces some errors when calculating the behavior of the system for times t <^ 
with a restriction of achievable cooling rates from above. Nevertheless, collective cooling can be much 
more efhcient than previously considered mechanisms j3- 

■^Note that the coupling constant g depends on the geometry of the respective setup. In case of a ring 
cavity, g and the annihilation operator c incorporate all possible modes of the quantised electromagnetic 
field in the resonator. 

"^We also observe that the cavity is in resonance with the atomic transition of Figure 1(a) and therefore 
the cavity does not couple to the vibrational modes of the particles. 



4 



non-resonant coupling of the laser Hamiltonian to the 1-2 transition of the particles. This 
term is negligible compared to the driving of the resonant excitation of the sideband with 
coupling strength y/N^riQ,^, if 0.,^ -C We neglect this non-resonant laser driving here since 
we do not expect it to have an effect on the conclusions drawn in the paper. It can only lead 
to an additional evolution between the states \l)p but cannot cause unwanted population 
outside the Dicke-symmetric subspace. Using the notation 

x, = ^/Nlr]n,, x = {J24y^^ = l^V^, b = J2ix,/x)K, y = VNg, (10) 

1/ V 

Eq. becomes 

i/j = ;ia;S+6 + %^+c + H.c. , (11) 

where the phonon annihilation operator h obeys the familiar commutator relation [6, 6^] = 1. 

The leakage of photons through the cavity mirrors is accounted for by considering the 
master equation [22] 

p = [Hi, p] + K {cpc^ - \ c^cp - \ pc^c) , (12) 

where p is the density matrix of the combined state of all particles, the common vibrational 
mode and the cavity field. Assuming regime (0) and restricting ourselves to Dicke states 
with / -C A^, spontaneous emission from the particles is negligible. Moreover, all noise terms 
(like heating) with amplitudes small compared to x, y and k can be neglected. 

The concrete form of Eq. (fT^ suggests that the stationary state pss corresponds to a 
coherent state of the form |a)p|/3)v|7)c with 5'~|q;)p = a\a)p, b\P)^ = and c\'y)c = 

7 |7)c. Indeed, the only solution of pss = is the state |0)p|0)v|0)c with all particles in the 
ground state, no photons in the cavity and no phonons in the common vibrational mode. 
Once the atoms have been initialised, the system loses its phonons within the time it takes 
to reach the stationary state. Since this time evolution is solely governed by the frequencies 
X = ^^/Nrjil, y = \/N g and k, we expect that this happens in the large limit in a time 
given by the smallest of these frequencies. 

To calculate the cooling rate explicitly, we derive a set of differential equations for the 
variables 

m = {b^b)p, n={c^c)p, S3 = (S'g)^ , (13) 

where m is the mean number of phonons with respect to the phonon number operator 6^6, 
n is the mean number of photons in the cavity mode and S3 relates to the mean number of 
particles in the excited state |1). It is also useful to consider the coherence quantities 

h = {S+b-S'b^)p, 
h = (5+c-5-ct),, 

h = {bc^ + b^c)p. (14) 
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The introduction of ()14|) is motivated by the fact that the system first builds up the co- 
herences ki, which, once established, provide an effective coupling between the different 
subsystems. Using Eq. (fT^ we obtain 

m = ixki , h = iyk2 — nn , S3 = — i {xki + yk2) . (15) 

The annihilation and creation of phonons and photons is accompanied by changes of the 
excitation number of the single particles. For / <^ A^, these fluctuations remains negligible 
on average and the approximations {S3b^b)p = Ssin, {S3C^c)p = s^n and {Ssi^bc^ c)) p = s^k^ 
can be adopted. We also neglect contributions of order one compared to N such as {S^S~)p. 
Then 

h = -j^i'^xm + yk3)s3, 

h = ~{2yn + xk3)s3-\nk2, 

ks = i{yki + xk2) - ^Kks . (16) 
These non-linear differential equations imply 

X 

m = — (k/c3 + 2^3) -{Kn + h). (17) 

2y y^ 

We see below, that the presence of a negative k^ and a positive n provides a cooling channel 
and plays a crucial role in the cooling process. 

Here we are interested in the cooling of a large number of particles. This allows us to 
solve the time evolution considering first the regime where k ^ and Eq. (jl7p becomes the 
conservation law 

2 

m--ks + %n = 0. (18) 

y y 

In the parameter regime (0) and given that m, n and S3 are of order N, the system reaches 
a stationary state with m, n and k^ constant on a time scale of the order 1/^/N. After 
such a time we might safely assume m = -h = ki = k2 = k^ = 0, and obtain from Eq. (fT^ 
and for k ^ that ki = k2 = 0, k^ = —{2x/y)m and n = {x'^/y'^)m. This stationary 
state is actually reached in a time given by the smallest among y/Ng and ^VNtjO., while 
K ~ controls longer-lived processes such as the one described by Eq. (fTTj) . Under the above 
assumption for the zeroth order in k, we find m = — [x^(x^ + y'^)/y^] to the first order in 
K. From this we get 

mit) = ttlq exp ( Kt] , (19) 

where itlq is the initial number of phonons in the system with respect to the above defined 
operator 6^6. The exponential decrease of the phonon population (see Figure |21^a)) amounts 
to the overall system cooling with a rate of the same order of magnitude as n. The result 
(fTIJ|) . which, we stress, holds under the condition of large A^, shows that the cooling of the 
system does not depend on the specific value of A^, and thus, provided A^ is large, it holds 
even in the case that not all the particles are initially prepared in their ground state, which 
helps the feasibility of the scheme. We also remark that such a behavior becomes possible 
only because of dissipation, namely the leakage of photons through the cavity mirrors. 



6 




Figure 2: Cooling evolution of common vibrational modes obtained from a numerical solution 
of Eqs. (fT3j) - (fTBj) (solid line) in comparison to Eq. (fTIJ|) (dashed line) for g = 10~^ k, rjQ = 
5 ■ 10~^ K, N = 10^ and mo = 10'^ (a). Cooling of individual phonon modes obtained from 
a numerical solution of Eq. (j25j) assuming n = k2 = and S3 = —^N for g = 10~^ n, 
r]Vt = b- 10"^ K,N = 10^ and mo = 10^ (b). 



4 Collective cooling of the individual motion of the 
atoms 

The collective regime established above may also be obtained in the case where each particle 
i couples to its own set of individual phonon modes 6,^ j. As above, we assume that the Rabi 
frequencies Vli, of the corresponding laser fields with detuning v are for all particles the same. 
Then the Hamiltonian of the system equals in the interaction picture 

Hi = Y, l^V^u (T+K,, + hg a+c + H.c. (20) 

u,i 

With the notation 

x, = VnIv^., x = iJ24y^^^ b^ = J2(^u/x)K,i, y = VNg, (21) 

where the hi obey the relation = Eq. (^01) simplifies to 

= ^ E ^ ^t^i + y ^tc + H.c. (22) 
ViV i 

Suppose that the system is initially prepared in a state with all particles in the ground state, 
the mean phonon number is the same for all particles and of about the same size as mo 
considered before and there are no photons in the cavity mode. Then the operator Y^f=i erf hi 
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has a similar effect on the system state as the operator S~^b in the previous case The 
net result is shown to be again collective cooling in the large N limit. 

Proceeding as before, we first calculate the stationary state p^s- Leakage of photons 
through the cavity mirrors is accounted for by using Eq. (fT^ . As before, the form of the 
Eqs. (fT^ and suggests that pss is the state with all particles in the ground state and no 
phonons and no photons in the cavity. Indeed, it obeys pss = 0. 

To calculate the cooling rate explicitly, we consider the expectation values 

m = N{b{h)p, h={c^c),, h = N{a3i)p (23) 

and the coherences 



k2 = VN (cr+c - (JiC^)p , 

h = N{bic^ + b\c)p, (24) 

where we sum over all particles. Using Eq. (fT^ with the Hamiltonian ()22j) and the same 
approximations as before, we obtain 

rh = ixki , h = iyk2 — ku , S3 = — i {xki + yk2) (25) 

and 

ki = --^ {2xm + yk3)s3 , 

h = {2Nyh + xk3)s3- ^Kk2, 

fcs = i{yki + xk2) - Inks (26) 

implying 

2 

fh = -^ {nh + 2^3) -^{nn + n). (27) 
2y y^ 

The differential equations and reveal that, for k ^ 0, the system reaches a stationary 
state within a time of the order one with /ci = ^2 = and k^ = —{2x/y) m, given m is of order 
N, while the cavity accumulates a small population of photons. However, n = {x'^/Ny'^) m 
remains small and the coherence k^ provides the main decay channel for the phonons in the 
system. From Eq. (P7j) we obtain for large N, in analogy to Eq. (fT^ . 

/ x'^ \ 

mit) = niQ exp ( Kt) , (28) 

V / 

where fho is the initial total phonon number. Thus, similar to the result (fT^ . the rate for the 
cooling of individual phonon modes is of the same order of magnitude as the cavity decay 
rate k, after a transition time given by the smallest among g and ^r]Q (see Figure El^b)). 
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5 Conclusions 



In this paper we propose a new cooling mechanism based on the collective excitation and 
de-excitation of particles trapped inside an optical cavity. In Section we showed that a 
large atomic sample in the presence of highy symmetric interactions, that treat all particles 
in the same way, behaves like a collection of bosonic particles. In Sections 0] and 01 we then 
analysed the two extreme cases, where the applied laser fields aim either at the cooling of 
only common or individual motions of the particles. For both cases, we predict cooling rates 
of the same order of magnitude as the cavity decay rate k. In the general case, one might 
argue that the number of vibrational modes is extremely large. However, we believe that 
similar cooling rates would be achievable in this case as well, as long as the initial phonon 
energy in the setup does not increase with the number of particles in the setup, which is 
in general not the case. 

We conclude with a few comments. The conservation law ()17|) obtained for k = 
allows to define the conserved quantity Q = Q with Q = m + (x^/?/^)n — {xjy) k^. The 
meaning of the time independence of Q is that one can shift the quantities ks, m and n by 
some constants without changing the dynamics of the system apart from the changes in the 
process of redistributing phonons, governed by fci, and their population m. The leakage of 
photons through the cavity mirrors [k ^ 0) disturbs the equilibrium expressed by Q = 
inducing a dynamical response, i.e. a quantum phase transition: The overall effect is the 
exponential decrease in the phonon population, namely cooling 

The treatment in this paper should be compared with the traditional approach where 
the attention is focused on the single particle behavior. In such a single particle treatment 
[21] the cooling depends exponentially on the cavity decay rate k. In this has to be 

about the same size as g and \'!]^^ which imposes a strong constraint on the time scale of 
the cooling process. The advantage of our treatment lies, on the contrary, in the fact that 
the numbers of particles introduces, in principle, the crucial freedom to tune k to be as 
large as \fNg and ^VNrjQ thus cooling the particles very efficiently. Although corrections 
might be necessary in realistic situations, a high cooling rate also reduces the number of 
collisions between particles, which certainly are a problem in less effective cooling schemes. 
We believe that considering the cavity-mediated collective field theoretical behavior opens a 
new perspective in particle cooling. 
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